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The phase transition of the quantum spin-f/2 frustrated Heisenberg antiferrofer- 
romagnet on an anisotropic square lattice is studied by using a variational treatment. 
The model is described by the Heisenberg Hamiltonian with two antiferromagnetic 
interactions: nearest-neighbor (NN) with different coupling strengths J\ and J[ along 
x and y directions competing with a next-nearest-neighbor coupling J2 (NNN). The 
ground state phase diagram in the (A, a) space, where A = J[/Ji and a = J2/J1, is 
obtained. Depending on the values of A and a, we obtain three different states: anti- 
ferromagnetic (AF), collinear antiferromagnetic (CAF) and quantum paramagnetic 
(QP). For an intermediate region Ai < A < f we observe a QP state between the 
ordered AF and CAF phases, which disappears for A above some critical value Ai. 
The boundaries between these ordered phases merge at the quantum critical endpoint 
(QCE). Below this QCE there is again a direct first-order transition between the 
AF and CAF phases, with a behavior approximately described by the classical line 
a c ~ A/2. 

PACS numbers: 75.fO.Jm, 05.30.-d, 75.40.-s, 75.40.Cx 



I. INTRODUCTION 



The study of the phase transition of frustrated spin systems on two-dimensional (2d) 
lattices is a central problem in modern condensed mater physics. A competition of ex- 
change interaction can lead to frustration, where spatial arrangement of magnetic ions in 
a crystal for which a simultaneous antiparallel ordering of all interacting spin is impossi- 



2 



ble. In particular, one of the frustrated 2d models most discussed is the quantum spin-1/2 
Heisenberg antiferromagnet on a square lattice with competing nearest-neighbor (NN) and 
next-nearest-neighbor (NNN) antiferromagnetic exchange interactions (known as J\ — J 2 
modeipffl. 

The criticality of this J\ — J 2 Heisenberg model on a square lattice are relatively well 
known at T = 0. There are two magnetically long-range ordered phases at small and 
at large values of a = J2/J1 separated by an intermediate quantum paramagnetic phase 
without magnetic long-range order in the region between a lc ~ 0.4 and a 2c — 0.6, where 
the properties of these disordered phase are still under intensive debate. For a < a\ c , the 
system possesses antiferromagnetic (AF) long-range order with wave vector Q = (tt, tt), with 
a staggered magnetization smaller than the saturated value (quantum fluctuations), which 
vanished continuously when a —> a.\ c . For a > a 2c we have two degenerate collinear states 
which are the helical states with pitch vectors Q = (tt, 0) and (0,7r). These two collinear 
states are characterized by a parallel spin orientation of nearest neighbors in vertical (or 
horizontal) direction and an antiparallel spin orientation of nearest neighbors in horizontal 
(or vertical) direction, and therefore exhibit Neel order within the initial sublattice A and B. 
At a = a 2c , the magnetization jumps from a nonzero to a zero value. The phase transition 
from Neel to the quantum paramagnetic state is second order, whereas the transition from 
the collinear to the quantum paramagnetic state is first order^LS Isaev, et alW^ have shown 
that the intermediate quantum paramagnetic is a (singlet) plaquette crystal, and the ground 
and first excited states are separated by a finite gap. 

The interest to study the two-dimensional J\ — J2 Heisenberg antiferromagnet have been 
greatly stimulated by its experimental realization in vanadium phosphates compounds^^S], 
such as Li 2 VOSi04, Li 2 VOGe04, and VOM0O4, which might be described by this frus- 
trated model in the case of J 2 — J\ (a = J%/ J\ — 1). These isostructural compounds are 
characterized by a layered structure containing V 4+ (S = 1/2) ions. The structure of V 4+ 
layer suggest that the superexchange is similar. In these compounds a second order phase 
transition to a long-range ordered magnetic phase has been observed. NMR spin-lattice 
relaxation measurements^ below T c shows that the order is collinear. Due to the two-fold 
degeneracy of the ground-state for a > 0.7 it is not possible to say a priori which will be 
the magnetic wave vector (i.e., Q = (tt, 0) and (0,tt)) below T c . On the other hand, such 
a scenario can change by considering spin-lattice coupling which will lift the degeneracy 
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of the ground-state and will lower its energy^. Then, any structural distortion should in- 
evitably reduce this competing interactions and thus reduces the frustration. In the case 
of this frustrated magnetic materials, the competing interactions are inequivalent but their 
topology and magnitudes can be tuned so that the strong quantum fluctuations destroy 
the long-range ordering. Experimentally the ground state phase diagram of frustrated com- 
pounds, described by the J\ — J 2 model, can be explored continuously from high to the low 
a = J2/J1 regime by applying high pressures (P), which modify the bonding lengths and 
angles. Recent results from x-ray diffraction measurements^ on the Li 2 VOSi04 compound 
has shown that the ratio a decreases by about 40% when the pressure increases from to 
7.6GPa. 

A generalization of the J\ — J 2 Heisenberg antiferromagnetic model on a square lattice 
was introduced by Nersesyan and TsvelilpS and studied by other groups^^l, the so-called 
Ji — J[ — J 2 model. In the Ji — J[ — J 2 model is considered inequivalence nn couplings Ji and 
J[ = A J\ in the two orthogonal spatial lattice dimensions with all the NNN bonds across the 
diagonals to have the same strength J 2 . Study of extensive band structure calculations^ for 
the vanadium phosphates ABVO(P04) 2 (AB=Pb 2 , SrZn, BaZn, and BaCd) have indicated 
four inequivalent exchange couplings: J\ and J[ between NN and J 2 and J 2 between NNN. 
For example, in SrZnVO(P04) 2 was estimated J[/J\ — 0.7 and J 2 /J 2 — 0.4 causing a 
distortion of the spin lattice. This spatial anisotropy tends to narrow the critical region and 
destroys it completely at a certain value of the interchain parameter A. 

On the other hand, by using the continuum limit of the J\ — J[ — J 2 spin-1/2 model 
Starykh and Balents 2 - 1 ' have shown that this transition splits into two, with the presence of 
an intermediate quantum paramagnetic (columnar dimer) phase for A < 1. Bishop, et aP^, 
by using coupled cluster treatment found the surprising and novel result that there exists a 
quantum triple point (QTP) with coordinates at (a t = 0.33 ± 0.02, \ t = 0.60 ± 0.03), below 
which there is a second-order phase transition between the AF and CAF phases while 
above this QTP are these two ordered phases separated by the intermediate magnetically 
disordered phase (VBS or RVB). The order parameters of both the AF and CAF phases 
vanish continuously both below and above the QTP, which is typical of second-order phase 
transition. There is some evidence that the transition between the CAF and intermediate 
phases is of first-order. Using exact diagonalizatiorP^ with small lattice of iV < 36 (6 x 6) 
size, the intermediate QP phase for all interval of A G [0, 1] has been obtained for the pure 
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spin- 1/2 Ji — J 2 model on a square lattice. These results are in accordance with results 
obtained by Starykh and Balentes^, that predicted not the QTP in the ground-state phase 
diagram recently observed by Bishop, et aZ™. 

The ground state (GS) properties of the two-dimensional frustrated Heisenberg antifer- 
romagnet have been investigated by various methods. The exact diagonalization starts from 
singlet states on pairs of sites, which cover the whole 2d lattice. However, the manifold of 
these states which can be constructed is nonorthogonal and overcomplete. This numerical 
methods are limited to small clusters iV < 6 x 6 = 36 due to storage problems. The com- 
putation on the largest cluster 6x6 has been performed by Schulz and co-workes^ 20 years 
ago. In spite of the great improvements achieved during this time, it is not possible so far 
to repeat this calculation for the next interesting cluster 8x8. This is only possible with 
other technique, as the quantum Monte Carlo simulation. Due to the progress in computer 
hardware and the increased efficiency in programing, very recently^ the GS of the quantum 
spin- 1/2 J\ — J2 model have been calculated by the Lanczos algoritm for a square lattice 
with N = 40 sites. 

The theoretical treatment of the frustrated quantum models is far from being trivial. 
Many of the standard many-body methods, such as quantum Monte Carlo techniques, may 
fail or become computationally infeasible to implement if frustration is present due to the 
minus-sign problem. Hence, there is considerable interest in any method that can deal with 
frustrated spin systems. This considerable qualitative difference in the ground state phase 
diagram in the a — A plane of the quantum spin- 1/2 J\ — J[ — J 2 model further motivates 
us to study this issue by alternative methods. 

Using a variational approximation, in which plaquettes of four spins are treated exactly, 
Oliveirapl has studied the ground state phase diagram of the pure J\ — J2 Heisenberg an- 
tiferromagnet on a square lattice, where the quantitative results are in good accordance 
with a more sophisticated method (exact diagonalization). In this work, we generalize this 
variational method to treat the anisotropic square lattice ( J\ — J[ — J 2 model). The rest 
of this paper is organized as follows: In Sec. II, the model is presented and a brief discus- 
sion of results. In Sec. Ill, the method is applied for the case of one plaquette with four 
spins interacting with other plaquette type mean field approximation. Main results will be 
presented in Sec. IV, as well as some discussions. Finally, in Sec. V we will give a brief 
summary. 
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II. MODEL 

The critical behavior of the quantum spin-1/2 J\ — J 2 Heisenberg model has been studied 
for many years, but very little has been done in the anisotropic square lattice case, which is 
described by following Hamiltonian: 

H — ^2 (Jl a i,j ' + J'l a i,j ' a i,j+l) +^2 H " Ci+l.j+l + &i+l,j ' j (1) 

(id) ((id)) 

where er^— (of*, cr^j, cr*J) is the spin-1/2 Pauli spin operators, the index labels the x 

(row) and y (column) components of the lattice sites. The first sum runs over all NN and the 
second sum runs over all NNN pairs. We denote the Hamiltonian (1) by J\ — J[ — J 2 model, 
with strength J\ along the row direction, J[ = XJi along the column direction, J 2 = aJ\ 
along the diagonals, and we assume all couplings to be positive with J[ < J\. 

The classical (S = 00) model (1) has only two ordered ground-states: AF (or Neel) for 
a > A/2 and columnar stripe (CAF) for a < A/2 separated by a first-order line at a c = A/2. 
Quantum fluctuations play a significant role in the magnetic phase diagram of the system 
at zero temperature. We will investigate the role of quantum fluctuations on the stability of 
the Neel and collinear phases. In the S = 1/2 case (quantum limit), the line splits into two 
phase transitions, where the ordered states (AF and CAF) are separated by an intermediate 
quantum paramagnetic (QP) phase, both on a square lattice. Exact diagonalizatiori-^ has 
estimated a critical line at cxcaf = oc c + for the transition between the CAF and QP 
states, and at oaf = ot c — ^ between the AF and QP states. The phase diagram in the 
a — A plane obtained is in accordance with Starykh and Balents^. However, the existence 
of QTP [quantum triple point) that was predicted by Bishop, et a/.^, is not present in 
their obtained phase diagram. Moreover, they found only presence of second-order phase 
transitions in the phase diagram. This contradictory qualitative results (existence or not of 
QTP) is the primary motivation behind this present work. 

On the other hand, a critical endpoint (CE) is a point in the phase diagram where a 
critical line meets and is truncated by a first-order line. This CE appear in the phase 
diagram of many physical systems such as binary fluid mixtures, superfluids, binary alloys, 
liquid crystals, certain ferromagnets, etc, and have been known for over a century^. Despite 
the CE long history, new singularities at the CE were predicted. Fisher and UptorP^ argued 
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that a new singularity in the curvature of the first-order phase transition line should arise at 
a CE. This prediction was confirmed by Fisher and Barbosa's^phenomenological studies for 
an exactly solvable spherical model. In conclusion of the analysis of the multicritical behavior 
observed in the ground-state phase diagram in the A — a plane for the J\ — J[ — J2 model, we 
have the presence of a quantum critical endpoint (QCE) and not QTP as mentioned other 
works^EH Therefore, the objective of this work is to obtain the QCE using the variational 
method, that was developed previously by OliveirapJ in the pure limit (A = 1) case. 



III. METHOD 



We first express the fluctuations around the classical ground state (AF and CAF phases), 

N/4 

where consider a trial vector state |^o) — Ti |0oz) f° r the ground state as a product of 

2=1 

plaquette state |</>oz) . We denote the plaquettes by I label, that is composed of four spins, 
where it do not overlap (mean field) on the square lattice as illustrated in figure 1. Each 
plaquette state is given by 



> = XX|7i>,, ( 2 ) 



l^oz 

n=l 

where {|1> = |+ - +-), \2) = |- + -+),|3) = |+ + — ),|4) = |- + +-),|5) = 

4 

I h+) , |6) — \-\ h)} is the vector basis with & z = Yl a t = 0' i a i} are re& l var i- 

i=l 

6 

ational parameters obeying the normalization condition Yl a n = !■ With this choice of 

n=l 

vector states, the mean value of the spin operator in each site of the plaquette is given by 
(<Ji) = rriiZ, rrii = (erf) = (0 O z \&f\ <poi), where the components in the x and y directions are 
null. 

Using the trial vector state defined in the Eq. (2), we obtain the magnetizations at each 
site that are given by 



mi = 2 (xu + yv — zw) , 
m 2 = 2 (—xu + yv + zw) , 
m 3 = 2 (xu — yv + zw) , 



(3) 
(4) 
(5) 



FIG. 1: Two-dimensional square lattice with a plaquette structure to be considered in this paper. 
The negrit plaquette is composed with the four 01, 02, 03, an d 04 spin operators that are considered 
inEq. (2). 

and 

m 4 = 2 {—xu — yv + zw) , (6) 

where we have used the same set of parameters (canonical transformation) of RefP^, i.e., x = 
(ai + a 2 ) / a/2, y = (a 3 + a 5 ) / \/2, 2: = (a 4 + a 6 ) / a/2, u = (ai - a 2 ) / a/2, u = (a 3 - 05) / a/2, 
and w = (a 4 — a^) / \/2, which obeys the normalization condition x 2 +?/ 2 +2: 2 +-u 2 +i; 2 +u> 2 = 1. 
The ground state energy per spin and unit of J±, E = (\5 ^0) /NJ\, is given by 

= -^01 + -^02, (7) 

with 



£ 01 = (^1.^2) A + A (^2-^3)^ + (^3-^4> A + A (^4^1) A + 

A (^> A . (^ 5 } B + A (J 2 ) A . tf 6 ) B + (b* 2 ) A . (^ 8 ) D + (^ 3 > A • (^ 9 > D , (8) 
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and 



E 02 = a [fafa) A + (^2.^ 4) A + fa) A ■ (#6) B + (^2) A ■ (^5)b + fa) A ■ fa) C + 

fa) A ■ (^9> D + fa) A ■ fa) D + (^6> fl • , (9) 

where (O) = (<f> 0ll \0\ 0o M ) is the mean value of a given observable (9 calculated in the 
vector state of the fi(— A,B,C,D) plaquette as illustrated in Fig. I. 

The variational energy can be evaluated using the properties of the spin- 1/2 Pauli oper- 
ator components, i.e., a z |±) = ± |±),<7 X |±) = |=f ) and a y |±) = ±i |=f), that is expressed 
for 



£0 ( ^1 + m 2 ) - 2 (A + 1) sV + S— ^ 



(y 2 + ^ 2 )-(, 2 + u; 2 ) 2 ] + 



2 (1 - A) (y 2 v 2 - z 2 co 2 ) + 2x (yX + z) + 



1 



— (y — z) z — v 2 — co 2 + Qx 2 u 2 — Qy 2 v 2 — Qz 2 uj 2 

(10) 



To obtain the minimum energy with a boundary condition given by normalization x 2 + 
y 2 + z 2 + u 2 + v 2 + w 2 = 1, we use the Lagrange multiplier method which correspond the 
minimization of the functional 



F{x, y, z, u, v, w,rj) — E — rj (x 2 + y 2 + z 2 + u 2 + v 2 + w 2 - l) . 



(11) 



The stationary solutions (5 T = 0) are obtained by solving the set of nonlinear equations 

- (A + 1) x - 4 (A + 1) xu 2 + 2 (yX + z) + 12axw 2 = 2r]x 
(1 - A) (y + 4yv 2 ) + 2xX - 2a (y - z) - \2ayv 2 = 2rjy 
- (1 - A) z - 4 (1 - A) zoj 2 + 2x + 2a (y - z) - \2azuo 2 = 2r]z 

- (A + l)u = 2r]u 
(1 — A) v — 2av = 2rjv 
— (1 — A) u — 2au = 2r)ui 

where 77 is the Lagrange multiplier. 



(12) 



IV. RESULTS 

The variational parameters x, y, z, u, v, w, and rj are determined simultaneously solving 
the system of equations (12) combined with the normalization condition x 2 + y 2 + z 2 + 
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u 2 + v 2 + w 2 = 1 for each phase. In the quantum paramagnetic (QP) phase we have 
mi = m 2 = m 3 = m 4 = 0. We note that in the isotropic limit (A = 1), our results reduce 
the same expression obtained by OliveiraP^. In this disordered phase, the ground state vector 
\4>qi) is an eigenvector of S^, where S/ is the total spin of the Zth plaquette of four spins, with 
zero eigenvalue (singlet state). In the AF ordered phase we have the boundary condition 
mi = —1712 = wi3 = — wi4, and in the CAF phase mi = vri2 = —m^ = —m^. 

The order parameters ttiaf — { m i~ m2 + m 3 — m 4 )/4 and tticaf = (mi + m 2 -wi3-m4)/4 
are numerically obtained as a function of frustration parameter a for a given value of spatial 
anisotropy A. We observe that the order parameter ttiaf goes smoothly to zero when the 
frustration parameter (a) increases from zero to atj. c (A) with A > Ai characterizing a second- 
order phase transition. A simple fitting of the form itiaf — («ic — otY m the vicinity of 
the second-order transition gives the same classical value for the critical exponent (3 = 1/2. 
On the other hand, for A > Ai and a > a 2c (A) the staggered magnetization m C AF increases 
monotonically with the frustration parameter a in the CAF phase, with a discontinuity of 
m CAF at a = Q!2c(A), which is a first-order phase transition. For A < Ai, the QP intermediate 
phase between the two ordered states (AF and CAF) disappears, and a direct transition 
between the magnetically ordered AF and CAF located at the crossing point a c ~ A/2 
correspondent to the classical value. 

The ground state (T = 0) phase diagram in the A — a plane is displayed in Fig. 2. 
The solid line indicate the critical points and the dashed lines represent first-order frontiers. 
We observe three different phases, namely: AF (antiferromagnetic), CAF (collinear anti- 
ferromagnetic) and QP (quantum paramagnetic). The AF and QP phases are separated 
by a second-order transition line «i c (A), while the QP and CAF phases are separated by 
a first-order transition line «2c(A). The presence of the interchain parameter A has the 
general effect of suppressing the QP phase. The QP region decreases gradually with the 
decrease of the A parameter, and it disappears completely at the quantum critical endpoint 
QCE=(Ai,«i) where the boundaries between these phases emerge. Below this QCE, i.e., 
for A < Ai, there is a direct first-order phase transition between the AF and CAF phases, 
with a transition point a c ~ A/2 (classical value). 

In order to illustrate the nature of the phase transition, we also show, in inset Fig. 2, the 
behavior of the staggered magnetization (order parameter) as a function of the frustration 
parameter (a) for A = 0.4(< Ai) and 0.80 (> Ai). From curves such as those shown in 
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Fig. 2 we see that for A = 0.8 there exists an intermediate region between the critical point 
a = ai c (X) at which m,AF(a) — > for the AF phase, characterizing a second-order transition, 
and the point a = «2 C (A) & t which the mcAF{&) order parameter presents a discontinuity 
for the CAF phase, characterizing a first-order transition. For A = 0.4, the order parameter 
of the AF phase decreases monotonically with increase of the frustration parameter from 
~ 0.81, for a = 0, to zero for a ~ 0.2 (~ A/2). In the CAF phase mcAF{.ot) decreases 
from Pd 0.95 for a = 1.0 to ~ 0.82 for a = a^d^) — A/2, characterizing a direct first-order 
transition between the magnetically ordered AF and CAF phases located at the crossing 
point. We note that the definition of the order parameter m M = (erf) = (0 O z |of | (j) i) (/j, = 
AF, CAF) difer of 1/2 factor when compared with calculations which use other methods (i.e., 
m M = (jS!f ) = (0oz |of I 0o«) /2). Therefore, in the limit of the not frustrated (a = 0) square 
lattice (A = 1) antiferromagnetic, solving the equations (12) and applying the corrections 
factor we found m^F — 0.41 which is consistent with the numerical results obtained by 
various methods such as series expansion, quantum Monte Carlo simulation, and others^, 
and can also be compared with experimental results for the K 2 NiF 4 , K 2 MnF 4 , and Rb2MnF 4 
comp ounds^H^. 



V. CONCLUSION 



In summary, we have studied the effects of quantum fluctuations due to spatial (A) and 

frustration (a) parameter in the quantum spin-1/2 J\ — J[ J 2 Heisenberg model. Using 

a variational method we calculated the sublattice magnetization for the AF and CAF 
phases. For values of A > 0.51 the frustration contributes significantly to the existence of a 
disordered intermediate state (QP) between the two AF and CAF ordered phases, while for 
A < 0.51, we have a direct first-order transition between the AF and CAF phases. We have 
observed, by analyzing the order parameters of the AF and CAF phases, that the phase 
transitions are of second and first-order between the AF-QP and CAF-QP, respectively. 
The obtained phase diagram can be compared with recent results which used effective-field 
theory^l anc i coupled-cluster method^, showing the same qualitative results predicting a 
paramagnetic region for small interlayer parameter (i.e., A > Ai), and for A < Ai this QP 
phase disappears by presenting a direct first-order transition between the AF and CAF 
phases. On the other hand, recent calculations of second order spin wave theory^ have 
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FIG. 2: Ground state phase diagram in the A — a plane for the quantum spin-1/2 J\ — J[ — 
Ji model on a square lattice, where a = Jij J\ and A = J[/J\. The dashed and solid lines 
corresponds the first- and second-order transitions lines, respectively. The black point represents 
the quantum critical endpoint (QCE). The notations indicated by AF, CAF and QP corresponds 
the antiferromagnetic, collinear antiferromagnetic and quantum paramagnetic phases, respectively. 
The dotted line correspond the classical solution a c = A/2. 

indicated that the intermediate QP phase exists for all A < 1 in accordance with results 
of exact diagonalizatiorM We speculate that by using a more sophisticated method, for 
example, quantum Monte Carlo simulations^ and density matrix renormalization group 
(DMRG) methocP^, this disordered region should disappear for certain values of A < A x . 
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